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Abstract 

We discuss the problem for detecting long-range correlations in se- 
quences of values obtained by generators of pseudo-random numbers. 
The basic idea is that the Holder exponent for a sufficiently long se- 
quence of uncorrelated random numbers has the fixed value h = 1/2. 
The presence of long-range correlations leads to deviation from this 
value. We calculate the Holder exponent by the method of multifractal 
detrended fluctuation analysis (MFDFA). We discuss frequently used 
tests for randomness, finite sample properties of the MFDFA, and the 
conditions for a correct application of the method. We observe that 
the fluctuation function F q used in the MFDFA reacts to trends caused 
by low periodicity presented in the pseudo-random number generator. 
In order to select appropriate generators from the numerous programs 
we propose a test for the ensemble properties of the generated pseudo- 
random sequences with respect to their robustness against presence of 
long-range correlations, and a selection rule which orders the genera- 
tors that pass the test. Selecting generators that successfully pass the 
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ensemble test and have good performance with respect to the selec- 
tion rule is not enough. For the selected generator we have to choose 
appropriate pseudo-random sequences for the length of the sequence 
required by the solved problem. This choice is based on the close- 
ness of the Holder exponent of the generated sequence to its value 1/2 
characteristic for the case of absence of correlations. 



1 Pseudorandom sequences, Holder exponent 
and multifractal detrended fluctuation anal- 



Computer random number generators have many applications in physics as 
for an example in the Monte Carlo methods |2| or in the nonlinear time 
series analysis [HJH|. Computers implement deterministic algorithms, hence 
they are not able to generate sequences of truly random numbers. Thus the 
randomness of the sequences is relative: What can be random enough for 
one application may not be random enough for another application. The 
computer generated sequences are called pseudorandom, and in order to be 
appropriate for scientific applications they have to satisfy many requirements 
and must pass extensive statistical tests. For an example the period of the 
generator (the number of different results before the generator repeats itself) 
must be as large as possible. Another requirement is that the generated num- 
bers should not be correlated among themselves. A significant class of such 
correlations are the long-range correlations which decay much slower than 
exponentially with time or distance and are observed in many systems in the 
Nature [SlIE]- The autocorrelation function c n = (r}iT)i +n ) of a pseudoran- 
dom sequence {r] a } ideally should be a Kronecker <5„ >0 - c n is not a convenient 
tool for detecting long-range correlations among the numbers in the sequence 
since the presence or absence of weak correlation at large n is usually masked 
by statistical fluctuations. Hence we have to apply a separate test to detect 
them. In this paper we propose such a test which is based on the simple fact 
that the Holder exponent for the sequence of uncorrelated random numbers 
has the fixed value of 1/2 (for more details see the Appendix). The deviation 
from this value for large n is evidence for possible long-range correlations in 
the sequences generated by the used pseudo-random numbers generator. 

The Holder exponent is a measure of irregularity of a curve at a singular 
point [HE!- Let us consider a function z(t). Its r-oscillation at value t is 
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The graph of z(t) is fractal if osc r (t)/r — > +00 uniformly with respect to 
t. For differentiable z(t) the ratio osc(i)/(2r) tends to dz/dt for r — ► 0. 
When the limes superior of the ratio is infinite, there is no derivative, and 
the Holder exponent measures the singularity of the graph of z(t) in this 
point. A function z(t) is Holderian of exponent h if a constant c exists such 
that for all t' 

(2) \z{t)-z{t') \<c\t-t'\ h 

or in terms of r— oscillations osc T (t) < cr h . If osc T (t) > cr h the function z{t) 
is ant i- Holderian of exponent h at the point t. The relation between c n and 
the Holder exponent for the case of Hurst noise is discussed in the Appendix. 
We recall from there the relationship (for large n) 

(3) c n = h(2h - l)n 2h ~ 2 . 

For the case of pure white noise (c n = 0), h = 1/2. When h > 1/2, c n > 0, 
i.e., the series of data are correlated and when h < 1/2, c n < 0, i.e., the data 
series are anti-correlated. 

We shall calculate the Holder exponent for our pseudorandom sequences 
by means of the MFDFA method [Oj which is briefly described in the Ap- 
pendix. We perform all calculations by means of the MFDFA(l),i.e., linear 
fit of the trend in each segment (When the fit is made by polynomial of p-th 
order the method is denoted as MFDFA (p)). Our choice p = 1 respects the 
fact that no trend is to be expected in the investigated sequences of numbers. 
Below we shall discuss a test for long-range correlations and not for a trend 
in the generated sequences of pseudorandom numbers. If such a trend exists 
it has to be detected by other tests. But even in the case when these addi- 
tional tests are not performed the MFDFA reacts to the presence of trends 
and nonstationarities. This sensitivity is know to exist for DFA |TTT|. ITT] and 
as we shall observe below it exists for the case of MFDFA too. 

We use as random number generators several generators from ^2] namely 
ranO, ran2, ran3, a quick and dirty generator which we call qdg as well as the 
generator based on the program G05CAF from the NAG library. ranO is the 
minimal standard generator which has to be satisfactory for most applica- 
tions but if its parameters are not appropriately chosen the generator can 
have correlations. The generator ran2 is claimed to be very good one and 
without serial correlations up to limits of its floating point precision. ran3 is 
based on the algorithm due to Knuth [T3]. qdg is based on a cycle containing 
the lines of code 
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jran=mod(jran*ia+ic,ini) 
ran=float(jran) /float (im) 

where im, ia, ic are appropriately chosen parameters. The above generators 
are not tested and compared with respect of long-range correlations. 

The paper is organized as follows. In the following section we discuss 
the properties of the investigated generators with respect for standard tests 
for randomness. Finite sample properties of the MFDFA method are in- 
vestigated in section III in order to select interval of appropriate values for 
the parameter q of the method for which the test for randomness has to be 
performed. In the section IV the behavior of the fluctuation function of the 
MFDFA method is discussed for the investigated generators and it is shown 
that the fluctuation function is a good indicator for deviation from fractality 
and presence of trends in the generators. In section V we discuss an en- 
semble test for long-range correlations in the random number generators as 
well as rules for selection of (i) most appropriate generator for the required 
length of the sequence and (ii) most appropriate sequences generated by a 
selected generator. Some concluding remarks are summarized in the last 
section where we discuss MFDFA(O) and MFDFA(l) with respect to their 
sensitivity to detect deviations from the case h — 1/2. In the Appendix we 
describe shortly the relationship between the correlation function and the 
Holder exponent as well as the MFDFA method. 

2 Behavior of generators with respect to stan- 
dard tests 

In order to be sure that ranO, ran2, ran3, qdg and G05CAF produce se- 
quences close to random ones we have tested them by standard tests such as 
a frequency test or a serial correlations test pQ. The general difficulty with 
statistical test on finite sequences lies in the statistical fluctuations. For the 
latter tests they are perfectly known. We expect that the generators gener- 
ate sequences of pseudorandom numbers with histograms uniform in some 
interval ( [0, 1] for ranO, ran2, ran3, G05CAF and [—1, 1] for qdg). In order to 
perform the frequency test we sort a generated sequence of N numbers into B 
bins with expected mean value of the numbers in each bin M = N/B. If the 
actual number of random numbers in the j— th bin is Mj we can construct 
the quantity 

(4) xl = (IMX 2 = O-MijjiMj - Mf/M}} 

3=1 
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with v — B — 1 degrees of freedom, xl must be different from because of 
presence of some fluctuations in a histogram of a finite sequence of numbers. 
However the value of xt must be not too large because large values are the 
evidence for the concentration of numbers in some bins and thus the generator 
is not random, x 2 was calculated for v = 49 degrees of freedom and for a 
good random number generator x 2 has the value of about 40 for a sequence 
of 50000 pseudorandom numbers. Characteristic values for the investigated 
generators are presented in Table I. We note (i) the good performance of the 
quick and dirty generator qdg for small length of the generated pseudorandom 
sequences opposite to the case of large length of the sequence and (ii) the fact 
that the different generators have different performance for different lengths 
of the sequence. We can conclude that all generators except the qdg passed 
this test and we observe that the choice of the appropriate generator depends 
on the lengths of the pseudorandom sequence we need. 

As a second test we have calculated the autocorrelation at lag v 

(5) c„ = — ((x n - (x))(x n - v - (x))) 

where (x) and a 2 are the mean and the variance of the sequence. Several 
typical results for the investigated generators are presented in Table II. The 
small values of the autocorrelations show that the generators successfully pass 
this test. This fact supports the observation that the long-range correlations 
could be quite good masked so that the pseudorandom number generators 
with such correlations could pass the standard correlation tests. We can 
conclude that except for the cases of very inappropriate generators the stan- 
dard tests do not supply us with much information about the question if 
the pseudorandom generator we want to use is free from long-range correla- 
tions. We have to design other tests to detect such correlations. The test 
and rules which are discussed below are based on the use of the Holder expo- 
nent which is calculated by means of the multifractal detrended fluctuation 
analysis (MFDFA) method. 

3 Finite sample properties of MFDFA 

Before we apply the MFDFA method to our pseudorandom number sequences 
we have to understand its finite sample properties. For infinite uncorrelated 
sequences of values the Holder exponent has the fixed value h — 1/2. We dis- 
cuss here in detail finite sample properties of sequences generated by two of 
the generators: the generator ran2 which is claimed to be quite good in [T2*] . 
and the generator ranO. The investigation of sequences generated by C05CAF 
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and ran3 leads to the same conclusions. We estimate the Holder exponent 
h(q) from a sample of length N for ensemble of pseudorandom sequences (For 
the role of the parameter q see the Appendix ). The MFDFA method will 
be consistent when h(q) — > 1/2 with increasing length N of the sample, and 
unbiased when the ensemble average (htm(q)) — > 1/2. Our observation is 
that the MFDFA method is consistent but biased for large | q |. For the in- 
vestigation of the finite sample properties of the two generators we have used 
at least 10 ensembles each of 25 random sequences (i.e. at least 250 different 
sequences). For larger number of sequences in the ensembles the results for 
h(q) do not change significantly and the only significant effect is the decreas- 
ing of the standard deviation of the mean. We test the consistency of the 
method by keeping parameters unchanged except for the length of the time 
series which is increased . Several characteristic results for the generator ranO 
are presented in Fig. 1. These results as well the results from other similar 
calculations starting from different seeds of the generator parameter idum of 
ran2 show that the MFDFA method is consistent. However, the method is 
biased for large | q | as it can be seen from Fig. 1 and Fig. 2. We observe that 
for small | q | (the numerical experiments lead us to the conclusion that small 
is | q | < 2) h(q) stays around 0.5 but for larger | q | significant deviations 
exist (i) when we increase the length of the sequence and keep s ma and s max 
constant (Fig.l), and (ii) when we keep the length of the sequence constant 
and increase s max (Fig. 2). A part of the bias of the method for large | q \ 
can be removed by increasing sample size because for small sample size we 
have limited number of values with low probability. Nevertheless our advice 
is that all test for long-range correlations must be performed for small | q \. 
As generators generate only pseudorandom sequences we can not expect to 
obtain exactly 1/2 as a value for h. Our numerical investigation has shown 
that differences of 0.005 are admissible, i.e., the generator (from the point 
of view of ensembles of generated sequences) can be considered safe with 
respect to long range correlations when for small | q | (between and 2) its 
ensemble h(q) spectrum is between 0.495 and 0.505. 

4 Behavior of the fluctuation function 

When we investigate sequences from a pseudo-random numbers generator 
first of all we have to check if the fluctuation function F q (s) from the MFDFA 
method really scales as a power of s for different values of q. We expect 
F q (s) to be a straight line on the log-log scale (panels (a),(b),(c)) of Fig. 3) 
. Deviations from this behavior are evidence for a problem. In panel (d) we 
see a broken line with a break point approximately at 1/2 of the period of 
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generator which is characteristic for sequences with periodic trend. Hence 
the MFDFA method reacts on trends in the generated sequences i.e. it can 
be used as a warning message for presence of trends too. In addition, even if 
the fluctuation function is a straight line, but the value of h is significantly 
different from 0.5 for small | q |, this is an evidence for presence of long-range 
correlations in the sequence of pseudo-random numbers. In panels (a)-(c) of 
Fig. 3 we observe that with increasing length of the sequence for ran2 F q (s) 
comes closer to a straight line. This could be expected because the larger 
length leads to a better statistics as the form of the sequence histogram 
approaches the ideal assumed histogram form. Panel (f) shows a comparison 
between fluctuation functions for characteristic sequences generated by ranO 
and ran2. We observe that ran2 is more robust than ranO with respect to 
long-range correlations because the fluctuation function lines for different q 
are closer to straight lines for the case of the generator ran2 for all range 
of segment lengths s. Fluctuation functions for ranO are more dispersed for 
large lengths of the segment s. 

We haven't obtained satisfactory results for the generator qdg. This 
means that the fluctuation function does not scale as power law of the length 
of the segment (as in the case of the vales of im=6075, ia=106, ic=1283), or 
the scaling exponent is quite significant from 0.5 as for an example in the 
case of the values of parameters im=233280, ia=9301, ic=49297. 

5 A test and selection rules 

On the basis of obtained results we can formulate a test and rules of selection 
among pseudo-random numbers generators and among sequences generated 
by a generator. The test determines the conditions under which a random 
number generator can be considered to generate large amount of sequences 
free from long-range correlations. We note here that no absolute test exists, 
i.e., a generator can be good for one required length of the sequence but 
another generator could be better for another length of the sequence. Let us 
have several pseudo-random numbers generators and let us want to choose 
these of them which generate large amount of sequences free of long-range 
correlations. This choice can be made on the basis of the following test on 
the ensembles of pseudorandom sequences. 

For a given random number generator take at least 10 different ensembles 
each of at least 25 pseudorandom sequences and calculate the Holder exponent 
by means of MFDFA(l) method for | q |< 2. If for all ensembles the fluctuation 
function F q (s) scales as power law for all q and the Holder exponent is between 
0.495 and 0.505 the random number generator can be considered to be able to 
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generate large amount of sequences free from long-range correlations. 

Generators which pass the above test have to be preferred with respect 
to generators that fail the test. Thus one can select generators with small 
probability of generating sequences possessing long-range correlations (for 
the required length of the sequence). The next question which arises is how 
to order the appropriate generators and to select one of them. The answer is 
given by a selection rule which is based on a power-law scaling property of 
the fluctuation function and on closeness of corresponding Holder exponent 
to 1/2 and states 

Let us have two pseudo-random numbers generators which pass the above 
test. Let us calculate the fluctuation function F q (s) for at least 10 different 
ensembles each of at least 25 sequences for the two generators. The generator 
which has closer to power law form of F q (s) for all q and for which the h(q) is 
closer to 1/2 is more robust with respect to long-range correlations. 

An extensive investigation leaded us to the following ranking of the gen- 
erators with respect of the test and the selection rule 

1. ran2 

2. G05CAF 

3. ran3 

4. ranO 

Despite the fact that a generator passes the test and it is selected by the 
selection rule each pseudorandom sequence has to be tested separately i.e. 
the generator can be chosen among a manifold of generators but nevertheless 
some of its sequences can have long-range correlations among their values. 
Therefore after choosing the most appropriate generator (for the required by 
the solved problem length of the pseudorandom sequence) we have to test 
every generated sequence for presence of long range correlations. The best 
sequences are selected by a selection rule which is analogous to the above 
selection rule namely that the appropriate sequences have Holder exponents 
most close to 1/2 for | q |< 2. Several results for pseudorandom sequences 
from different generators are presented in Table III. The sequences are chosen 
for illustration of the fact that each generator can generate sequences for 
which h is considerably different from 1/2. 

6 Concluding remarks 

In this paper we have used the multifractal detrended fluctuation analysis to 
investigate the behavior of the Holder exponent for sequences of pseudoran- 
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dom numbers obtained by several random number generators. Theoretically 
the Holder exponent for a large enough sequence of random numbers must 
have a fixed value h = 1/2 regardless of the order q of the fluctuation func- 
tion in the MFDFA method. The deviations from randomness lead to three 
kinds of changes: (i) The fluctuation function F q (s) is not a straight line on 
a log-log plot. This is evidence for presence of some trend in the generated 
sequences i.e. the generator is very bad one; (ii) When F q (s) scales as a 
power law of s, h(q) is a straight line significantly different from 1/2 for all 
values of q for | q |< 2. This indicates presence of long-range correlations but 
if h do not depend on q the generated sequence has monofractal properties 
up to smallest investigated length of the segments of the MFDFA; and (iii) 
h(q) could be close to 1/2 for some values of | q | but not for all values for | q | 
for | q |< 2. This means presence of long-range correlations and multifractal 
properties of the generated sequences. The existence of bias at large | q \ 
means that MFDFA must be used very carefully when one calculates charac- 
teristic fractal quantities for sequences of values with multifractal properties. 
Our case here is a monofractal one (theoretically h has a value, independent 
on q) but nevertheless we have to take into account this bias by restriction 
on the values of | q \. 

The simplest variants of the MFDFA from the point of view of the fitting 
polynomial are MFDFA(O) and MFDFA(l). Above we have used MFDFA(l). 
It is possible to use MFDFA(O) in the formulated test and selection rules (i.e. 
to use MFDFA without local detrending). Let us discuss the properties of 
MFDFA(O) and MFDFA(l) with respect to their application to sequences of 
pseudorandom numbers. Let us assume that our sequence of pseudorandom 
numbers is divided into N s segments each of length s and let us write the 
profile function Y for the z/-th segment as 

(6) Y[(u-l)s + i]=i^ 2 + 6 v (i) 

where S u (i) is the deviation of Y from i 1 ^ 2 at the z-th value of the u-th. 
segment. For the fitting polynomial we assume 

(7) y„(i) = a v % + b v 

where a v and b v are constant coefficients. If a u = b u = we have MFDFA(O) 
variant of the MFDFA method. When a v ^ and (or) b v ^ we have the 
MFDFA(l) variant of the MFDFA method. 

From @ and (|IJ) we obtain for the variation F 2 {y 1 s) for the i/-th segment 
(for large enough values of s) 

F 2 (u, a) = 8 -- a u b u s + - £ Sl{%) + - £ 2t 1/2 5 u (z) - 

1 S 2 = 1 S i = l 
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a u s 3 / 2 - Ks 1 ' 2 

(8) 

Let for simplicity below q — 2. Then for the fluctuation function we have 



(9) F 2 (s) = ^P + Q + S 

where 

1 2Ns fs 

do) p = _ S (__ m „ s 

contains the terms corresponding to the ideal case, 

(11) Q = Q* + Q** 
where 

i s 2N S 

(12) Q* = ^££K« + 2* 1/2 <^)] 

(is) Q« = — y: ep^w - 2 w)] 

Z6JV s i=l l/=l 

are the terms containing the fluctuations 5„(i). Q* depends only on the 
fluctuations and Q** depends on the fluctuations and on the trend. Finally 

contains the terms which depend only on the local trend. 

For the case of MFDFA(O) (a„ = K = 0) P = Q** = 0. in the ideal case 
Q* should be negligible and -^(2) oc s 1 / 2 for large s i.e. the Holder exponent 
is h — 1/2. In the real case Q* could lead to deviation from h = 1/2 and if 
these deviations are large this is an indicator for presence of problems in the 
generated sequence. 

For the case of MFDFA(l) the sensitivity can be higher as we have more 
terms that can affect the value of the Holder exponent. The term P in 
this case is again proportional to s 1 / 2 independent on the values of the local 
trend coefficients a u and b u . For large enough s and large enough sequence 
of random numbers all other terms should be negligible and h oc 1/2. In 
the real situation deviations can come from Q*, Q** and from R which for 
large s is proportional to s 2 . When the local trends are correlated (i.e. some 
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kind of global trend is presented) then R ( which for large s could become 
larger than P) could lead to h close to 1 instead to h ~ 1/2 as in the case 
without long-range correlations. Similar situation is observed for other values 
of q. For an example when q = 4 the term P is a sum of terms of the kind 
s 2 (l/4 + 2a 2 b 2 u + a u h u ) (i.e. P oc s 2 ) and the dominant for large s member 
of R is proportional to s 4 . In summary if we do not want to use MFDFA 
with local detrending we can base our test and rules on MFDFA(O). In this 
case the deviation from the ideal case (i.e. from h — 1/2) is evaluated on the 
basis of quantities like Q*. When we use MFDFA(l) we can gain additional 
sensitivity with respect to deviations of h from 1/2. 

Finally we note that the MFDFA is not the only possibility for calculation 
of Holder exponent. Another method is the WTMM (wavelet transform 
modulus maxima) method El El El E]- MFDFA is shown to have 
slight advantages for negative q and short time series For long time 
series the WTMM could have advantages with respect to MFDFA and thus 
for length of sequences larger than 10 7 values the formulated tests above 
should be used on the basis of the results of the WTMM method. In such a 
case instead of scaling of the fluctuation function F q (s) one has to study the 
scaling of the partition function Z q (a) used in the WTMM. 
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A The importance of the exponent h and MFDFA 

Let us consider the general one-dimensional random walk in discrete time. It 
is a sum of steps which can be either discrete or continuous. Let the walking 
particle starts at the origin. After N steps its position X N is a sum of N 
mutually independent random variables. We are interested in a situation 
when these variables have the same distribution function F of mean \x and 
finite variance a 2 . It can be shown on the basis of the central limit theorem 
that for large N Xn is approximately normally distributed with mean Nfi 
and variance Na 2 PUimini- Thus the standard deviat ion of this random 
walk scales as N h with h = 1/2. 

In order to understand better the Holder exponent h let us consider a 
sequence of observations r]i, i = 1, 2, , . . . , N. We choose a reference size m 
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and sizes m such that mp = N, where p is the number of segments each 
of size m. For fixed m we calculate the mean value M(m) and the stan- 
dard deviation S(m) for each segment and construct the relative dispersion 
R(m) = S(m)/M(m). From general manifold of all possible data sequences 
we shall consider these for which we can observe 



(15) [R(m)/R(m )] = {m/m 



,h-l 



where for simplicity we consider h to be a constant. Mandelbrot |7| calls 
similar fluctuations and noise Hurst noise in order to notify the important 
contribution of Hurst [22] to the research of processes with long-range corre- 
lations. By means of (fTB^l we obtain a system of equations for the correlation 
functions 

(16) C T = (ViVi+r) = C x jC 2 

where 

N-t N 

(17) G x = [N/(N - r)] £ VlVt+T - (£ Vif/N 

i=l i=l 
N N 

(is) tf 2 = 5>? - (5>) a /tf, 

1=1 i=l 

for r = 1, 2, . . . , n. The system is 



n-l 



:n 2h -n 



(19) E(^-l)^ = (l/2)( 

i=l 

and it has the solution [231 121] 

(20) c n = (l/2)n 2h 1(1 + (l/n)) 2h - 2 + (1 - (1/n)) 



For very large value of n we can represent the term in [. . .] as Taylor series 
and as a result we obtain the relationship between the autocorrelation and 
the Holder exponent 

(21) c n = h(2h - l)n 2h - 2 . 

Recently a multifractal detrended fluctuation analysis method has been 
proposed for the analysis of long-range correlation of nonstationary time 
series jSJ. Here we present the variant of the method useful for calculation of 
significantly different from zero positive Holder exponents. The first step of 
the method is to calculate the mean (x) of the investigated time series. Then 
we calculate the profile function = jy k=1 (xk — (x)), i = 1, 2, . . . , N. After 
calculation of Yi we divide the time series into segments and calculate the 
variation for each segment. The division is into N s =int (iV/s) segments and 
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because the obtained segments would not include some data at the end of the 
investigated time series, additional N s segments are added, which start from 
the last value of the sequence in the direction to the first value of sequence. 
In order to calculate the variation we have to calculate the local trend (the 
fitting polynomial y v {i) for each segment of length s where s is between an 
appropriate minimum and maximum value). The variations are defined as 

(22) F 2 (y, s) = - ± {Y[(u - l)s + i] - y u (z)} 2 

S 8=1 

for the first N s segments and 

(23) F\v, «) = -£; {Y [N -{y- N)s + i] - y^)} 2 

for the second N s segments. Finally we construct the q-th order fluctuation 
function 

2N S 

(24) F q (s) = {[l/(2N s )]^[F 2 (u,s)]^}^. 

For monofractal time series F q (s) has to scale as s of constant power h which 
for sequences of random numbers has the value 1/2. Even in presence of 
local correlations extending up to a characteristic range s* the exponent 
h = 1/2 would be unchanged when s » s*. If the correlations do not have 
characteristic lengths the exponent h would be different from 1/2 fIE\ I2fi] . 
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Table 1: \ 2 test for several pseudorandom number generators. idum=2 for 



ranO, ran2, ran3. im= 

Number of points 


=6075, 
ranO 


ia=106, 
ran2 


ic=1283 for qdg. 
ran3 qdg G05CAF 


10 4 


46.30 


39.83 


45.13 7.14 42.25 


10 5 


46.29 


44.38 


57.13 2.35 41.34 


10 6 


42.04 


43.29 


38.84 16.54 59.68 



10 7 " ~ 41.25 57.33 41.38 165.84 56.35 



Table 2: Autocorrelations c v for several pseudorandom number generators. 
Length of generated sequences is 10 5 values. idum=2 for ranO, ran2, ran3. 



im=6075, ia=106, ic=1283 for qdg. 
v ranO ran2 ran3 


qdg 


G05CAF 


10 2 -0.0061 -0.00098 -0.0013 


0.0017 


0.00086 


10 3 -0.0008 -0.0033 -0.0076 


-0.00009 


0.0018 


10 4 -0.0029 0.0043 -0.0061 


-0.000098 


0.0055 



5-10 4 -0.0045 0.0013 -0.0023 -0.0048 0.0019 



Table 3: Holder exponent h for several pseudorandom number generators. 
Length of the pseudorandom sequences: 10 6 values. *: initialized with 
G05CBF(0). **: initialized with G05CBF (2). 



generator 


q=-2.95 


q=-1.95 


q=-0.95 


q=0.05 


q=1.05 




ranO (idum= 


= 10) 


0.500 ±0.005 


0.501 ±0.004 


0.502 ±0.004 


0.503 ±0.004 


0.505 ±0.003 


0. 


ranO (idum= 


=20) 


0.543 ±0.004 


0.542 ±0.004 


0.540 ±0.004 


0.539 ±0.004 


0.537 ±0.004 


0. 


ran2 (idum= 


= 10) 


0.536 ±0.005 


0.535 ±0.005 


0.535 ±0.005 


0.535 ±0.005 


0.535 ±0.005 


0. 


ran2 (idum= 


=20) 


0.533 ±0.007 


0.534 ±0.006 


0.534 ±0.006 


0.535 ±0.005 


0.535 ±0.005 


0. 


ran2 (idum= 


=4) 


0.489 ±0.006 


0.489 ±0.006 


0.489 ±0.006 


0.489 ± 0.005 


0.489 ± 0.005 


0. 


ran3 (idum= 


= 10) 


0.513 ±0.005 


0.512 ±0.005 


0.510 ±0.004 


0.508 ±0.004 


0.506 ±0.004 


0. 


ran3 (idum= 


=20) 


0.470 ±0.004 


0.469 ±0.004 


0.469 ±0.004 


0.470 ± 0.004 


0.470 ± 0.004 


0. 


ran3 (idum= 


=2) 


0.499 ±0.005 


0.499 ±0.005 


0.500 ±0.005 


0.501 ±0.004 


0.502 ±0.004 


0. 


G05CAF 


* 


0.479 ± 0.007 


0.479 ± 0.006 


0.480 ± 0.006 


0.481 ±0.005 


0.482 ±0.005 


0. 


G05CAF ' 


<* 


0.494 ±0.005 


0.495 ± 0.004 


0.495 ± 0.004 


0.496 ±0.004 


0.497 ±0.004 


0. 
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Figure 1: Influence of the length of sequence on the h(q) spectrum. Generator 
ranO, ensemble averages for ensemble of 25 time series. The segment length s 
is between s m i n = 10 and s max = 1000 for all panels. The lengths of sequences 
are 10 4 points for panel (a) 10 5 points for panel (b) and 10 6 points for panel 
(c). The h(q) spectra are denoted by solid lines and the dashed lines denote 
the correspondent h(q) spectrum plus/minus the standard deviation of the 
mean for the correspondent value of q. 
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Figure 2: Influence of the length of segment s on the h(q) spectrum of the gen- 
erator ran2. The investigated sequences of numbers are obtained for idum=2 
and have length of 250000 values. The h(q) spectrum is denoted by a solid 
line. The two doted lines denote the h(q) spectrum plus/minus the standard 
deviation of the mean for the correspondent value of q. Panel (a): h(q) spec- 
trum for s between s min = 10 and s max = 1000. Panel (b): h(q) spectrum 
for s between s min = 100 and s max = 10000. Panel (c): h(q) spectrum for s 
between s m i n = 1000 and s max = 25000. 
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Figure 3: Fluctuation function F q (s) for the pseudorandom sequences from 
ran2 and ranO. From bottom to the top at each panel the lines are for q = 
2,4,6,8. Panels (a)-(c): Results for ran2, idum=2. Panel (a): Fluctuation 
functions for a sequence of 10 4 points. The values of the Holder exponent 
are: h = 0.505 ± 0.002 (q = 2); h = 0.505 ± 0.004 (q = 4); h = 0.503 ± 0.004 
(q = 6); h = 0.499 ± 0.005 (q = 8). Panel (b): Fluctuation function for 
a sequence of 10 5 points, h = 0.499 ± 0.001 (q = 2); h = 0.501 ± 0.001 
(q = 4); h = 0.505 ± 0.002 (q = 6); h = 0.509 ± 0.003 (q = 8). Panel (c): 
Fluctuation function for a sequence of 10 6 points, h = 0.500 ±0.0004 (q — 2); 
h = 0.502 ± 0.001 (q = 4); h = 0.504 ± 0.001 (q = 6); h = 0.506 ± 0.001 
(q — 8). Panel (d): Fluctuation function for ranO with parameters a = 343, 
q = 127773, r = 2836, m = aq + r = 43828975 . Period of the generator is 
4880, idum=2. Panel (e): Fluctuation function for random number generator 
ranO with parameters a = 16807, q = 127773, r = 2836, m = aq + r = 
2147483647. idum=2, sequence of 10 6 values, h = 0.511 ± 0.003 for q = 2; 
h = 0.506±0.003 for q = 4; h = 0.501 ±0.002 forg = 6;/i = 0.498±0.002 for 
q = S. a = 16807, q = 127773, r = 2836, m = aq + r = 2147483647 for ranO. 
Panel (f): Comparison between fluctuation functions for sequences generated 
by the generators ranO and ran2. Length of the sequences 10 5 points. idum=2 
for the two generators. Parameters of the generators a = 16807, q = 127773, 
r = 2836, m = aq + r = 2147483647 for ranO. Circles: Fluctuation functions 
for ran2. Squares: Fluctuation functions for ranO. Solid lines: Power-law fits 
for the sequence generated by ran2. From bottom to the top: q — 2, 4, 6, 8. 
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